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Abstract 

We study stress relaxation in a strongly segregated lamellar mesophase of diblock copolymers. We consider the 
extreme limit in which chains are highly stretched and with their junction points confined to narrow interfaces. A 
lamella can be divided into "stress blobs" at some distance z from the interface, with well defined local modulus 
G(z,uj) at frequency w. For sliding (compressional) stress the total modulus is transmitted in series (parallel) across 
the layer. We evaluate the local G{z,t) which shows, for a given height, a very broad spectrum of relaxation times. 
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I. INTRODUCTION 



Block copolymers made with incompatible blocks form a rich variety of domain structures. At equilibrium, a dense 
collection of monodisperse A-B diblock copolymer chains assemble themselves in minimum free energy configurations. 
Symmetric diblock copolymers with highly repulsive interactions between A-B segments microphase separate into a 
lamellar order with A-B junction points highly confined into interfaces to decrease the number of contacts between 
the two different blocks. The width of these lamellae is determined by a balance between the reduction of unfavorable 
contacts between A and B blocks and the decrease of chain entropy due to the stretching of the block polymer. 

The rheological behavior of diblock copolymers in a lamellar mesophase is of great importance due to the numerous 
practical applications of these systems. The understanding of how polymer domains flow and of how rheological 
properties differ from those of the much more familiar homogeneous melt presents an important and challenging 
theoretical problem. Exploratory rheological studies of lamellar domains have revealed strong departures form the 
behavior of simple melts. |^ Theories have begun to grapple with these observations. @ @ @ The motion of end- 
confined chains appears to differ strongly from that of simpler polymers. Thus, the dynamic response to deformation 
in block copolymer domains shows distinctive retardation and broadening; |^ correspondingly, the disentanglement 
dynamics of end-confined chains are expected to be much different from that of unconfined polymer melts and solutions. 
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Two cases for the linear response to a small step strain should be considered [[7| : a compression of the layers and a 
shear tending to slide the layers past one another. The block chains have their translational motion highly restricted 
as one chain end is strongly confined to the interface, so usual reptation, j|] which is believed to be the dominant 
relaxation mechanism of homogeneous system consisting of long ffexible polymers chains such as melts (and also, for 
example, the disordered phase of diblock copolymers) is highly suppressed. The disentanglent motion of the block 
chains with one end anchored at the interface resembles that of star polymers, Q for which relaxation proceeds by 
arm retraction along the tube (contour length fluctuations). In addition, the block chains are highly stretched along 
the direction perpendicular to the lamella, while still retaining their unperturbed (Gaussian) size in the directions 
parallel to the interface. In fact, relaxation proceeds at different rate in different places (depending how far the regions 
of interest are from the interface) and, for example, one can define a local stress modulus G{z, t) with z the distance 
from the interface, with regions at different z involving quite different time scales. In this contribution we investigate 
these distinctive and novel features of heterogeneous materials. 



II. THE EQUILIBRIUM LAMELLA 



Here we restate several relevant properties of copolymer chains as given, for example, in Refs. |^ and Q For 
simplicity, we consider that all interfaces are flat, parallel to each other and separated by a distance 2h. We suppose 
that each surface contains a chains per unit of area. In the incompressible melt state each chain fills a volume V 
proportional to the chain molecular weight (and to its chemical length). Hence, the thickness of the layer h is just 
given by ft, = Va. Polymers in the melt satisfy Gaussian statistics. Thus, the end-to-end distance, which is the 
characteristic dimension of an unperturbed polymer coil, is simply given by = V/h, with b a microscopic length 
which only depends on the local structure of the chain and liquid, usually of about a few Angstroms. We only study 
here the long chain limit of V >> 6^. 

The block chains ending at the interfaces fill uniformly the regions between planes so each chain must extend over 
distances of the order of h. Stretching the chain a distance h reduces the entropy of the random walk thus requiring 
work of the order of {h/RY ksT. We are interested in high energy interfaces with (h/R) >> 1. If a step strain is 
applied to the liquid, at very short times it behaves like a rubber with modulus Gq. From the elastic energy stored 
by the liquid one can obtain the typical volume per entanglement R^ via Go = ksT/R'^. This Re is independent of 
the molecular weight of the chain] hence long chains (i? >> Re) are highly entangled. 

Since the copolymers are strongly stretched in the lamellar mesophase, configurational integrals are dominated by 
the classical extremun of the action (Hamiltonian) , which corresponds to the most probable ("classical") 

configuration of the block extending from the interface into a microdomain. In this approximation, the free energy 
functional for the copolymer blocks in a mesophase domain consist of non-interacting Gaussian chains in a self- 
consistent external potential (pressure field) p{z), where z is the distance perpendicular to the interface. ||ll| [|| This 
p{z) is the free energy cost of bringing a segment to the height z from some initial reference position, and it forces 
the chains to assume a stretched state. We call v the volume displaced by the piece of the chain extending from the 
free end at zo to the height z. The "classical" configuration of a chain is given by the trajectory z{v) of a fictitious 



Newtonian particle of unit mass in an external potential p{z) which starting at some point zq at "time" v = falls to 
the z = surface (interface) in a "time" V. The "equal time" property for monodisperse systems gives a parabolic 
potential p{z) = po{l — (h'^/z^)) with po/ksT — n'^h'^b/SV^ and the chain "equation of motion" 
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It remains to specify the distribution of end positions zq. This distribution is found by imposing the constitutive 
properties of the melt or solution under consideration. For example, for a marginal solvent the local volume fraction 
is proportional to the pressure p. Instead, for a melt the local volume fraction is unity everywhere and we obtain 



e(zo) 



azo 



h^[l - (z2//i2)]l/2 



(2) 



where e(zQ) is the distribution of free ends per unit area per unit height. Chains from opposite interfaces can penetrate 
each other as long as work done against the potential p{z) is of the order of ksT. Although the energy p^V of these 
chains is typically much greater than ksT there is a small region near the midplane where p << pq. In a region of size 
^ — R{R/h)^'^/'^ near this midplane chain segments of size ^ have energy of order fc^T or less. In this "interpenetration 
zone" the effects of stretching are unimportant and the chains from the two opposite interfaces interpenetrate freely. 
In the large V limit of interest, we have h >> R >>£,>> Re, so segments that live in the interpenetration zone are 
still highly entangled. Figure 1 illustrates the relevant lengths of the equilibrium lamellar phase of strongly segregated 
diblock copolymers. 



Length scales controlling entanglement 

+ 




stress relaxes when 
most chain ends have 
passed through 



FIG. 1. Relevant lengths in a lamella. 



III. DISENTANGLEMENT PROCESS IN LAMELLAR DIBLOCK COPOLYMERS 



After a step-strain is applied to the liquid, the initial shear distorts the entangled copolymer chains. Since the 
number of entanglements per unit of volume is a local property of the melt, the initial value of the associated modulus 
Go is the same as that of a homopolymer melt. Over time the chains disentangle and the associated distortion 
dissapears. For the relaxation of sliding stress, the relevant entanglements are those of the interpenetration region. 
The stress is initially held by the entanglements in the interpenetration zone as chains from one interface exert force 
on chains on the other. Sinc;c there is only weak interpenetration, one would expect that the two opposite sides 
disengage relatively fast. On the other hand, over scales larger than Rg the stress redistributes more or less uniformly 
propagating down to the interfaces, so not only the interpenetration region is of relevance in the relaxation mechanism 
of the sliding stress. Compressional stress directly involves the whole chain and not only interpenetration segments. 
These entanglements are typically between adjacent chains of the same layer, so disentanglement of opposing layers 
is not sufficient to relax the stress. Thus, compressional relaxation is expected to procede slower than sliding stress 
relaxation. 

At the most local level, stress is transmitted along the chains, but in a larger scale, within a radius several times i?e, 
chains "collide" with many other chains redistributing the stress more of less uniformly among adjacent chains. At 
this distance the stress can be treated as a smoothly varying quantity. We then divide the liquid into into imaginary 
cubes ("stress blobs") each with a well defined dynamic modulus G{z,u!) at a height z and frequency u. To evaluate 
the local modulus G{z, uj) we imagine subjecting the single stress blob to a step-strain experiment illustrated in 
Figure 2. The stress relaxes as the initial entanglements disappear. Notice that while the stress is defined locally, the 
disentanglement motion is not. To release an entanglement between two chains, it is required that the end of either 
chain pass through the entanglement point. Of course, the ends are typically many stress blobs away. Knowing the 
chains trailing volumes v from the height z to their free ends provides one of the necesary ingredients to predict how 
quickly a given stress blob relaxes. 



FIG. 2. A blob of size of a few times Re is subject to a stress-strain experiment. Although the stress is defined locally, the 
free ends of the entangled chains could be far away. 

As ends of the block polymers are attached to the interface, reptation of the chains is highly restricted, and relaxation 
of stress points is expected to proceed by retraction of the tube in analogy with star polymers. The time distribution 




for the retraction of segments of volume v provides a second important ingredient to evaluate the local modulus 
G{z, t). Given the relaxation behavior of these stress blobs, one can readily find how the overall stress relaxes. Layers 
must be treated as a composite material of many sublayers, each with a dynamic modulus G{z,uj). For the sliding 
shear, the stress is transmitted in series across the layer. During the sliding stress relaxation, the stresses are constant 
for all z while the strain is the average of the strains over all z. The weakest element dominates (blobs closer to the 
height h). In the compressional stress, the system is qualitatively like a stack of rubber sheets which are clamped 
together and stretched: the stress is transmitted in parallel along each sheet. The strongest moduli dominates and 
the relaxation time is controlled by the blobs closest to the interface. Clearly the relaxation is strongly anisotropic. 
The local modulus G{z,t) is then given by 
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The quantity pz (v) is the (normalized) density of probability that a monomer in the stress blob at a height z belongs 
to a chain that, from that height to the free chain end, displaces a volume v. Pv{t) gives the probability that a trailing 
segment of volume v has not yet retraced its length in a time t. Then, the integral in Eq. (3) gives the survival 
probability at time < of a stress point at z, assuming that entanglements are two body processes in accord with the 
double reptation model In our simplest picture we also assume intra- and inter-blob statistical independence. 
That is, the release of a stress point in a blob at height z is independent of similar event at z', and that releasing a 
constraint within a stress blob will not trigger the release of many others. We also assume that entanglement points 
are uniformly distributed (as monomers in the melt are). 



A. Height Dependent Trailing Volume Distribution 



The distribution Pz{v) is obtained within the above sketched "classical approximation" of the self-consistent field 
method appropriate for weak excluded volume and moderately high surface coverage for strongly stretched chains. The 
local melt volume fraction 4>z{zo) at a height z contributed by chains with free ends in the interval Azo = [20, z^ + dzo] 
(with z < zq < h) is obtained as 

(f)z{zo)dzo = f'"'f] dzo , (4) 
I dz/dv I 

where e{zo)dz() given by (2) is the number of fictitious particles per unit area falling to the interface at z = from 

the height interval Azq and \ dz/dv \ is the particles velocity at z. 

For a given height z, when the free end is at zq, the volume v of the chain trailing segment is uniquely determined 
by the equation of motion, Eq (1), which relates z, zq and v. Then, the density of trailing volume at z is simply 
obtained from (4) changing variables from zq to v and keeping z fixed. That is, we should replace zq by zo(z, v) and 
multiply (4) by the Jacobian dzQ —\ dz^/dv \ dv. Then, the density distribution of trailing volume u at a height z is 
given by 

^y)dv = '-^^TTT- I dzo/dv I dv . (5) 
I dz/dv 

After performing the algebra and using equations (1), (2) and (5) we obtain 

Pz{v) = ''-^ TJ^ ,0<V< V,nax{z), (6) 
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with Pz{v) vanishing for Vmax{z) < v < V . The function Pz{v) posseses a integrable singularity at Vmax/V = 
{2 / ■n)arcos{z / h) indicating that most chains in the stress blob at height z have a trailing segments of volume Vmax{z). 
Clearly, near the middle of the lamella, when z w h, the trailing segments are very short. Conversely, near the 
interface with z ~ 0, the trailing segments are of the size of the full chain. Nevertheless, and in contrast with 
the Alexander [Q-de Gennes [|l^ picture, there is a distribution of trailing volumes which, in turn, gives rise to a 
distribution of relaxation times (or frequencies) over an extended time region. These results which arise from the 
"classical mechanical" description of the stretched polymer chain are slightly altered when "quantum" fluctuations are 
taken into consideration. (The major modification are the disappearence of the singularity in pz{v) and the existence 
of the penetration zone of width ^). 



B. Trailing Volume Survival Probability 



The calculation of the probability Pv{t) proceeds similarly to that of Pearson and Helfand in their study of 
relaxation of star polymers. In the tube model, the primitive chain consists of L/a segments with L = V/hRf, the 
length of the primitive tube and a = Re its step length (i?e is the typical distance between entanglements) . 

In the long chain limit under consideration there are a large number a of entanglements per chain, that is a = 
/ Rl ^ 1. We define a variable x = L — L which measures how much the primitive chain has contracted. Pv{t) is 
the probability that in time t the primitive tube of initial length L+AL, with AL the root-mean-square fluctuation 
of the tube length, has not yet contracted an amount y = {v/V)L by an spontaneous thermal fluctuation. AL — 
(L/(3a)^/^) « {tfiDy/"^ where D = ksT/Q, ( and tn. are, respectively, the diffusion constant, the chain friction 
constant, and the Rouse time. 

We call p{x, t; xq; y) the density of probability that the end of the chain is at position x at time i, given that initially 
{t = 0) was at I xo |< AL and provided that the length of the primitive tube was never smaller than L — y. It can be 
obtained by solving a flrst-passage problem | p6| . This probability density is the same as that of a diffusing particle 

(the free chain end) in the harmonic potential U — {3akBT/2L )x'^ with an absorbing barrier at y. Then, the survival 
probability at time t of a tail of volume v is given by Pv{t) = dx p{x, t; Xo;y), since, as long a.s x < y, the primitive 
chain has never contracted by an amount larger than y. In the limiting case oit^tfi and of x ^ AL the probability 
Pv{t) is obtained as 



(7) 



where the time ty is the typical time at which a segment of trailing volume v (measured from the free chain end) has 
retracted the tube, and it is given by 



tv = tv 
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The time ty at which the whole chain has retraced the tube is 
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(27r/3a)i/2 exp(3a/2). 



(8) 



(9) 



hence, it is exponentially large in the chain molecular weig ht M {a = R'^/Rl = M/Me, with Mg the entanglement 
molecular weight). Clearly, for long chains tv is much larger than both, the Rouse time tfi oc M^, and the entanglement 
time te oc M"^'^ of linear melts and concentrated solutions. Similar conclusions apply to ti, if a ^ 1 and v > V/a^^'^. 
Eq. (8) indicates that ty+Au is several times larger that ty if Av ~ l/a. This very strong dependence of ty on the 
volume V of the trailing segment suggest that the stress relaxation of the lamellar system is very sensitive to the 
trailing volume distribution. Moreover, from the above discussion one can conclude that, for large values of a, the 
fimction Py{t) is very well approximated by the theta function 



Py(t) = Q{v - v{t)) , when a > 1, 



(10) 



with v{t) obtained by inverting ty of (8). This is because, from Eq.(8), a given time t' has associated to it a volume 
v' such that ty> = t' . Then, if v > v' {v < v'), we have ty >> ty> {ty << ty>). 

Now we are ready to evaluate the local modulus G{z, t) defined in (3). Using Eqs. (6) and (10) we obtain 
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which is valid for tu << t < tv ■ This integral can be solved in closed form as 
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The logarithm of the local modulus \ogG{z,t) as a function of log(t/ty) is depicted in Figure 3 for various values of 
(z/h). Clearly, the trailing volume distribution leads to a very broad (unbounded in the large a limit) spectrum of 
relaxation times at a given height z above the interface. 




!og(t/tv) 

FIG. 3. log(G(z,t)/Go) vs. \og{t/tv) for various values of {z/h). We have taken (2a/3) = 100. 
IV. DISCUSSION AND CONCLUSIONS 



Stress relaxation in heterogeneous systems such as the lamellar phase of diblock copolymers proceeds at different 
rates in different spatial regions as shown in Figure 3. Thus we find that at intermediate times the diblock layer 
has three co-existing zones of strongly different rheology: the region near the midplane is a quickly relaxing liquid; 
the regions near each A-B interface are unrelaxed rubber. Remarkably, this heterogeneity arises from different parts 
within a given chain. We would have obtained similar heterogeneity even if we had neglected the dispersion of the 
chains through the layer, as in Alexander |lj]-de Gennes model. When we include this dispersion, as in Figure 3, 
we find in addition a very broad relaxation spectrum, even at a given height z. This broad spectrum results because 
of a) a distribution of tail lengths w at a given height z, and b) the exponential increase of the relaxation time with 
the tail length v. 

Contrary to our expectations, Eq.(12) predicts that G{z,t) = at the midplane (z = h). This is because within 
the "classical approximation" we have neglected the penetration zone of width ^. In this region, chains are Gaussian 
and hence the distribution of trailing segments is exponential in the trailing volume v; that is Pz=o{v) = exp(— w/w^), 
with = the mean volume of trailing segments living in the penetration region. In this region, segments 
are still liighly entangled » Re) and stress relaxation proceeds as discussed in Sec. 3 replacing the trailing 
volume distribution (6) by the exponential one. Then, we obtain the local modulus for a stress blob at the midplane 
G{z = 0,<) = {kBT/Rl)eK-p{-v{t)/v^) valid for tR«t< t,^. 

The G{z,t) of (12) is the mean value at time t of the stress modulus for a blob at height z. In fact, stress moduli 
at (z, t) are distributed according to the binomial distribution. This is because in a given blob at a given time, there 
are many stress points, each having a probability p to be alive. The probability p{z,t) is given by the factor in 
Eq.(ll) that multiplies the plateau modulus Go. In principle, the local modulus G(z, t) could be measured using local 
probes. On the other hand, the local modulus is the main ingredient to evaluate the total modulus for either sliding 
or compressional stress. It is straightforward to find these moduli using Eq. (12). It is also straightforward to extend 
this equation to account for the interpenetration zone. Our work in this direction is in progress. 
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